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We propose a general one-dimensional continuous formulation to analyze the vibrational modes of 
antenna-like nanomechanical resonators consisting of two symmetric arrays of cantilevers affixed to 
a central nano-beam. The cantilever arrays can have arbitrary density and length profile along the 
beam. We obtain the secular equation that allows for the determination of their frequency spectrum 
and illustrate the results on the particular examples of structures with constant or alternating 
cantilever length profiles. We show that our analytical results capture the vibration spectrum of such 
resonators and elucidate key relationships that could prove advantageous for experimental device 
performance. Furthermore, using a perturbative approach to treat the nonlinear and dissipative 
dynamics of driven structures, we analyze the anharmonic coupling between two specific widely 
spaced modes of the coupled-element device, with direct application to experiments. 
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I. INTRODUCTION 

Mechanical resonators, especially in the micro- and 
nano-range, are currently used in many research areas to 
investigate fundamental physics problems such as ultra- 
sensitive force and mass detection [1L single spin detec- 
tion , gravitational wave detection 3[ or quantum mea- 
surement and computation to cite a few. A key aspect 
of their performance is ultra-high resonance frequency of 
vibration - demonstrated up to several gigahertz as the 
nano-structures are made from stiffer materials such as 
UNC diamond or silicon carbide , and designed with 
submicron critical dimensions. The coupling of arrays of 
high frequency resonators can yield enhanced response 
characteristics in hybrid resonator designs, and as such 
holds primary interest in a plethora of technological ap- 
plications. Detailed analytical and numerical modeling of 
coupled element resonators is critical to the understand- 
ing and utilization of such devices. 

In a previous publication 7j we presented both dis- 
crete and continuum models for a class of coupled-element 
nanomechanical antenna-like structures that have been 
experimentally realized 0, H, Q . The continuum model 
was shown to capture the linear modal response spec- 
trum of the device in close agreement with finite element 
simulations. As a sequel to the previous work, here we 
generalize the analysis to coupled-element structures with 
arbitrary cantilever lengths and distributions along the 
central beam. By judicious choice of cantilever array dis- 
tribution, the response spectrum of the structure can be 
engineered for specific performance characteristics. As an 
illustration, we treat the particular cases of constant and 
alternating cantilever length arrays. Motivated by ex- 
perimental observations of modal coupling phenomena, 
we further expand the analysis to include weak nonlinear 
and dissipative contributions to the dynamics of a driven 
cantilever-beam and derive the amplitude-frequency re- 
sponse that describes the anharmonic coupling taking 



place between two of its specific widely spaced modes. 



II. ARBRITRARY CANTILEVER 
DISTRIBUTION 



We consider a doubly-clamped elastic beam with two 
symmetric arrays of cantilevers with lengths l(x) attached 
on both sides, see Fig. [T] We denote the beam out-of- 
plane deflection by y(x, i), x G [0, L] and by rj(x, £, t), the 
deflection with respect to y(x, t) of the cantilevers located 
at point x on the beam. Restricting our investigation 
to modes symmetric with respect to the beam, we have 

fe[o,i(4 

With subscripts b and c referring to quantities related 
to the beam and the cantilevers, respectively, we de- 
note the rigidities by £b and £ c , the masses per unit 
length by fib and fi c and the widths by Wb and w c . If 
we assume that all elements have the same thickness, 
£b/£c = Mft/Mc = Wb/w c [7j. The cantilever density is 
defined as p(x) = dN(x)/dx, where N(x) is the total 
number of cantilevers between and x. Assuming that 
shear deformation and rotary inertia of the structure are 
negligible, we work within the Euler-Bcrnouilli beam ap- 
proximation. Although not essential to our approach, 
this assumption greatly simplifies the interpretation of 
the results. More sophisticated beam theories can for in- 
stance be found in Ref. 11J and references therein. The 
Euler-Bernouilli Lagrangian of the structure depicted in 
Fig. Q}: is given by 
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FIG. 1: Schematic (2D) top views of mechanical antenna 
structures with symmetric arrays of cantilevers (along £-axis). 
a) Constant lengths, b) Alternating lengths, c) Length profile, 
l(x), and density along the central beam (z-axis) are arbitrary. 
Deflections occur perpendicular to the (x, £)-plane. 
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Cantilever boundary conditions apply to the deflection 
rj(x, £, t) and read 
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Boundary conditions for the beam can be arbitrarily 
specified. For a clamped-clamped beam, for instance, 
they would be y(0,t) = y'(0,t) = and y(L,t) = 
y'(L, t) = 0. From the least action principle, SS[y, rj\ = 0, 
where S[y,rj\ = JC[y,Tj](t) dt, we obtain the following 
equations of motion 



d 4 y(x,t) d 2 y(x,t) d 3 n(x,^t) 
Z-b — V 515 = -£cPK%) 
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The right hand side of Eq. is the shear force exerted 
by the cantilevers along the beam. Eq. ([¥]) combines 



the equations of motion of the cantilevers in a single one 
via T]{x, £, t). Its right hand side is the inertial driving 
force exerted by the central beam on the cantilevers. The 
energy of system ©-(jU), whose expression is obtained 
from Eq. (JTJ) by changing the sign of the elastic energy 
terms, is a constant of the motion. 

We now look for mode solutions of the form 

y(x, t) = Y(x) cos tut ; rj(x, £, t) = H(x, £) cos Lot. (5) 

Defining the useful quantities a and 7 by 

- 4 - '"' •- • - - al(x), (6) 
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we can solve Eq. Q together with boundary conditions 
© to obtain 

H(x,0 = [A+( 7 )cos(a0 +A+( 7 )sinK)+ 

A{ (7) cosh«) + A 2 - (7) sinh«) - 1] Y(x) (7) 



where 



Aid) 



1 + cos 7 cosh 7 =p sin 7 sinh 7 



± 



2(1 + cos7COsh7) 
cos 7 sinh 7 + sin 7 cosh 7 



(8) 
(9) 



2(1 + COS7COSI17) 

Using Eq. © again and Eqs. ®-©, we can evaluate 
the right hand side of the equation of motion for the 
beam, d 3 H/d£\= = -2a 3 A+ (~f)Y(x) and Eq. © can 
finally be cast into the form (the superscript (iv) denotes 
a fourth derivative) 

Y^ lv \x) - [V(a;x] + a 4 ] Y{x) = 0, (10) 
where the "potential" V(a;x) is given by 



V(a;x) 



w c a cos 7 sinh 7 + sm 7 cosh 7 

— p(x)a . 

Wb 1 + cos 7 cosh 7 



(11) 



The frequency spectrum of the antenna is obtained by 
solving Eq. (jTUJ) with four boundary conditions for the 
beam deflection Y(x) which results in a secular equa- 
tion for the parameter a that is, for the frequency u> 
(see Eq. ©). For example, if cantilevers and beam 
alike are assumed to be ID-like, the density reads p(x) = 
2 5(x—Xj), where the Xj's are the locations of the can- 
tilevers along the beam. The peculiarity of Eq. (jTUJ) is 
to have an x-dependent potential whenever the density or 
the length profile vary along the beam. Unless an explicit 
solution is known, it can be solved by projection onto an 
orthonormal set of modes satisfying appropriate bound- 
ary conditions for the beam. Eq. (fit)]) then becomes an 
homogenous matrix equation whose determinant provides 
the secular equation that yields the frequency spectrum 
of the structure. 



III. CONSTANT AND ALTERNATING 
CANTILEVER LENGTH ARRAYS 

Constant cantilever length - We now turn to the spe- 
cial case depicted in Fig. QJ, where two dense arrays of 
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FIG. 2: Frequency spectra of a constant cantilever length 
structure (left) and of an alternating cantilever length struc- 
ture (right). Frequencies to n ^ are normalized to the funda- 
mental's, ui.i. Physically relevant levels (such that n < TV) 
are represented as solid lines and others as dashed lines. Other 
parameters are given in the text. 



N regularly spaced cantilevers with constant length I are 
located symmetrically with respect to the central beam. 
As previously stated, in a one-dimensional setting, their 
exact density is twice the sum of N equispaced delta func- 
tions. When N is large enough, it is physically relevant 
to replace the density by its average value over the beam, 
p = 2N/L, as long as the modal shapes of the beam under 
investigation wash out the discreteness of the cantilever 
distribution. This requires for instance that the number 
of nodes of the said modal shapes, n, be small compared 
to the number of cantilevers on one side, N. Hence the 
condition of validity of our approach, n <C N. 

The physical meaning of replacing the exact density by 
its average is to spread the local shear forces exerted by 
the cantilevers on the beam around their support in or- 
der to form a continuous force density. This results in a 
cantilever continuum represented by the continuous de- 
flection T)(x,£,t). Notice that, although bi-dimensional, 
the cantilever continuum is distinct from an elastic plate 
because adjacent cantilevers arc not elastically connected. 
As p(x) is now replaced by p and 7 = al is a constant, the 
"potential" V(a; x) does not depend on x anymore and 
the solution to Eq. (fTO)) is straightforward. Its modes are 
given by 



Y n (x) = Aip n (x/L), 



(12) 



where A is an overall amplitude and {(p n (u)J, an or- 
thonormalized set of modes satisfying 

1 

<Pn V H u ) = Pt,<Pn(u), [ duip n (u)ip m (u) = 8 m>n . (13) 


The Pn's are the positive roots of the secular equation 
obtained from the boundary conditions for the central 



beam. For instance, for a clamped-clamped (+) or a 
clamped-free (— ) beam, 



COS /3 n COSh (3 n = ±1. 



(14) 



The secular equation for a, conveniently expressed here 
in terms of 7 = al, is obtained from Eqs. (|10[) and (|T2|l 

as 

2Nm c ,, cos 7 sinh 7 + sin 7 cosh 7 , ( l[3 n 

7 ^7 

TOfc 1 + cos 7 cosh 7 \ L 

where nib = f-bL and m c = p c l denote the mass of the 
beam and of a cantilever, respectively. For each particu- 
lar value of p n , Eq. (|15jl has an infinite number of roots 
7n,fcj k > 1, from which the frequency is evaluated as 




(16) 



The corresponding modes are determined from Eqs. (|12|) 
and where 7 is replaced by its value y n ,k- As we see 
from Eq. (112p . the modal shape of the beam is the same 
for the antenna as it would be for a beam without can- 
tilevers attached, a result also confirmed by finite element 
analysis 7]. This surprising fact originates in the special 
form of the force density exerted by the cantilevers on the 
beam that is proportional to the beam deflection Y n (x), 
thereby producing no change in the modal shape even 
though the frequency spectrum is modified by the first 
term of Eq. (JT5J) for N > 0. 

Eq. (|15[) has three relevant parameters: the length 
ratio, A = l/L, the product of the cantilever number 
by the width ratio, v = 2Nw c /wb and fi n . It can be 
shown that, for all n > 1 and k > 2, its solutions satisfy 
7oo,fc-i < 7n,fc < 7oo,fc where 700^ satisfies the secular 
equation of a simple cantilever, 



cos7oo fc cosh 7 



00, & 



1 = 0. 



(17) 



Notice that although the values 700, fc are parameter in- 
dependent, the corresponding frequencies, w^^, derived 
from (Tl7j|) still depend on I and other physical quantities. 
The frequencies satisfy 0Joo.k-i < ^>n,k < Woa,ki which 
explains the generic "band" structure of the spectrum il- 
lustrated in Fig. [2] (right) on the specific example of an 
antenna clamped at both ends with N = 20 cantilevers, 
tu c = Wb and l/L = 0.05. 

To gain some insight in the fine structure of the bands, 
we observe that Eq. (|15[) can be solved approximately for 
large n and for small or large A. These limits all require 
its first term to become large for finite values of 7, which 
leads to j n j. = + A„^ fc with 
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A^ <oo,fc + th 
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(18) 



In this expression, t^k = tan7 00)fe , th^^ = tanh7 oc ^ fc ; 
= fc if A„. fc < and fc = fc + 1 if A„ ifc > 0. Eq. (JI5J| is 



ft, — n, 11 i_\„ jfc ^ u aim ft, — ft, ■ 

valid for |A„, fc | < 1, k > 1 1 
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FIG. 3: Fundamental, 71.1 (red), and first collective, 71,2 
(green), levels versus the length ratio A = l/L for N = 
5, 10, 20, 50, 100 and 500. For clarity, N = 5 and N = 500 
are displayed as dashed and dotted lines, respectively. 



We can show that regardless of the boundary condi- 
tions, A nj fe vanishes as n becomes large and the frequen- 
cies w n ,fc cluster at cjoo.fc, the upper edge of band k. Pro- 
vided A is small enough, Eq. (fT8|) indicates that an accu- 
mulation of frequencies also occurs at the lower edge of 
band k + 1, a phenomenon observed in Fig. [5] In this 
limit, the frequency band-gap, Auik = <^i,fc+i — ^>oo.k is 
given by AuJk — 2w o,fcAi i fc/7oo,fc- It is seen to become 
constant as k increases. Clearly, from Eq. (fTH|) . if n is 
fixed and k > 2, 7„.fc tends to 7oo,fc-i (7oo,fc) when the 
length ratio A tends to zero (infinity). This indicates that, 
for band 2 and higher, the cantilever dynamics (see Eq. 
([T7|) ) dominate the motion of the structure in these two 
limits. This is not true however for the first (fundamen- 
tal) band whose lower edge tends to zero as A — > and 
whose dynamics is governed by the beam's in this limit, 
a mass loading effect. 

An interesting example of level motion is given in Fig. 
[3]where the fundamental mode, 71,1, and the first collec- 
tive mode, 71,2, are plotted against the length ratio A for 
various cantilever numbers, N. The central beam is again 
clamped at both ends and the width ratio is w c /wb = 0.5 
(v = N). For all N, 71.1 interpolates between zero and 
7oo,i as A increases and we note that, for a given length 
ratio, the fundamental frequency always shifts downwards 
when cantilevers are added to the central beam. In the 
small A regime, this is the manifestation of a mass loading 
effect while for large A, this is predicted by Eq. (fT8|) . Sim- 
ilarly, 71.2 interpolates beween 700,1 and 700,2 but there is 
a special ratio, A^ 2 , f° r which it becomes N -independent 
(curve intersection in Fig. [3]). For A < A\ 2 , 71,2 increases 
with N and decreases otherwise, as predicted by (IT8l) for 
small and large A, respectively. 

Such a ratio exists in fact for any level "f n ,k, k > 2. It is 
given by A* k — (i 2k _ s / {2f3 n ) where is solution to the 
clamped-clamped beam equation (fl"4")) while j3 n is solu- 



tion to the actual beam secular equation. For that ratio, 
7n,fe = ( l Pn)/L. Using S b /fJ, b = S c / p c and Eq. HI]), the 
corresponding frequency is thus w* k = \J Zbj Pb(Pn/ L) 2 , 
which is precisely the nth frequency of a beam without 
cantilevers attached! For the clamped-clamped antenna 
of Fig. El \* 12 = 1/2 and uj* 12 = ^/ £ b / Vb{Pi / L) 2 . Thus, 
if the cantilever length is half the beam's, the 1st collec- 
tive frequency of the structure does not depend on the 
number of cantilevers attached and is equal to the fun- 
damental frequency of the central beam alone. Although 
this conclusion is drawn from our continuum approach, it 
has been verified for the dynamics of the exact 3D struc- 
ture solved by the finite-element method. 

We also note that, all other parameters being fixed, the 
density of states (number of states per unit frequency) 
increases with N and that, for k > 2, the lower band 
edges generically shift upwards. The opposite occurs 
however for the first band (fundamental group) whose 
lower edge shifts downwards as N increases. The role 
of the cantilevers is indeed different for the fundamental 
group (k = 1), where they essentially produce a mass 
loading of the central beam, than for higher bands where 
their dynamics take over and shift the antenna frequen- 
cies upwards. 

Alternating cantilever length arrays - We briefly com- 
ment here on the case of arrays of cantilevers with al- 
ternating lengths, l\ and l 2 , h < h (Fig. QJ)). Using the 
continuum approach, we obtain the corresponding secular 
equation 

7 4 + ^ 2 ^ (7) + ^ 2 ^ (e7)= fh^Y (19) 

where 7 = al%, e — I2/I1 and where a and A 2 arc defined 
by Eqs. j6|) and (|9|), respectively. w\ and W2 refer to the 
widths of the cantilevers and p\ — 2Ni/L, p 2 — 2N 2 /L to 
their average density, Ni being the number of cantilevers 
with length /, on one side of the structure. Notice that 
for h =l 2 , Eq. (HI]) reduces to Eq. (JTSJ) as it should. The 
spectrum calculated from Eq. (fig)) is displayed in Fig. [5] 
(left). For comparison with its constant length counter- 
part, we have assumed w\ — u>2 = wt, Pi — P2 = N/L 
with N = 20 and h/L = 0.05 while e = 0.8. The spec- 
trum of the alternating structure shows additional bands. 
They are due to the term A 2 (£7) that becomes infinite 
when 7 = 7oo,fc/e 7 bands ensuing as previously discussed. 
This example shows how the frequency spectrum of a 
more complex antenna structure can be easily understood 
and provides a way of engineering frequency spectra with 
specific features useful in experiments. 

IV. NONLINEAR SYSTEM DRIVEN BY TWO 
FREQUENCIES 

A. Motivation 

It has been observed experimentally [u| that, if in ad- 
dition to being driven at the frequency of the fundamental 
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mode the antenna is also driven at the frequency of the 
collective mode, the frequency peak of the fundamental 
mode experiences a slight shift. This frequency shift is 
the signature of a modal coupling that occurs because of 
the presence of nonlinearity and dissipation in the sys- 
tem [l3[ . To explain the interaction of these two widely 
spaced modes, we investigate the effect of a two- frequency 
driving on the response of a weakly nonlinear and dissi- 
pative antenna structure. We supplement the equations 
of motion of our continuum model, Eqs. ([3]) and (UJ), with 
nonlinear terms that take into account the possible ma- 
terial and geometrical nonlinearities of the structure and 
with damping terms proportional to the transverse veloc- 
ity of the elements. These terms, small compared to the 
amplitude of vibration of the antenna, are treated within 
the Lagrangian approach described in [13| to derive the 
frequency-amplitude relation of the model. 



B. Lagrangian approach 

It is convenient at this stage to define the mode (n, k) of 
the constant cantilever length structure by the 2-colunm 
vector ip nik = (y n (x, t), r) n ,k(x, £, t)) T . It involves two pa- 
rameters, p n and 7n,fc, respective solutions to the secular 
equations Eqs. (fT4|) and (fT5|) . The corresponding fre- 
quency u> n ,k is determined by Eq. (fTo| while its modal 
shape is given by Eqs. ([7]) and (fL2"| . In particular, the 
fundamental mode is given by the parameters (/3i,7i,i) 
and the collective mode we are interested in by the pa- 
rameters , 71 2) • For clarity in the notations, we re- 
name these two modes tpi = (yi(x,t),rju(x,£ i ,t)) T and 
'02 = (U2(x,t),rj2i(x,£ : ,t)) T and denote their respective 
frequencies by u>i and u> 2 where 



Pf) 2 , *6{1,2}. 



(20) 



Experimental data (to be published elsewhere) shows 
that ipx and ip 2 are coupled in the sense that the 
amplitude-frequency curve (resonance) of the fundamen- 
tal mode is modified when the higher mode is driven. 
This coupling is attributed to the presence of nonlinear- 
ities in the resonator structure. To explain this phe- 
nomenon, we treat the nonlinearities and the damping 
affecting the system as a perturbation of the fundamen- 
tal and collective modes, ipi and tp 2 . These perturbative 
terms are responsible for a modulation of the linear modes 
that we evaluate by a multiple scale method. Hereafter, 
we closely follow the Lagrangian approach of Ref . [HI be- 
cause it offers a particularly suitable framework to derive 
the modulation equations. 

When the driving amplitude (or power) is small 
enough, the solution ip — (y(x,t),r](x,£ l ,t)) T to the non- 
linear equations of motion is, to a good approximation, 
given by a superposition of ip\ and ip 2 with slowly modu- 
lated amplitudes. Following the multiple scale approach, 



we write ip as 



+A 2 (T 2 )i Pl (u)e 



iu! 2 T 



+ ex. 



(21) 



In this expression, c.c. denotes the complex conjugate, e 
is a small bookkeeping parameter, u = x/L, v = 
Xk{v) — Hi t k(x,^)/Yi(x) and <p\{u) is defined in Eq. 
(flU|) . According to the multiple scale method, we have 
introduced two time scales, Tq = t and T 2 = e 2 t. The 
complex slowly varying modulation amplitudes are de- 
noted by Ai(T 2 ) and A 2 {T 2 ). 

The Lagrangian of our system is given by 



£ 



dx 







2\dt 



2 [dx 2 







2 I dt dt 



2 



d 2 r\ 



(NLT) + (Fi cos(Oi) + F 2 cos(wt)) / ydx, (22) 



where (NLT) stands for "Nonlinear Terms" . To express 
the fact that the driving frequencies, f2 and u>, are close 
to the linear frequencies of the fundamental and excited 
modes, we write them as 



CI = u>\ + e 2 CTi ; oj = uo 2 + e 2 a 2 . 



(23) 



To describe the nonlinear response of the cantilevers and 
the central beam, neglecting the effects of rotary inertia 
and shear deformations, we add the following nonlinear 
terms to the Lagrangian [l3[ 



NLT= / dx 





dy 
dx' 



dx' 



£b fdy_ d 2 yY 
2 \dx dx 2 J 



1/(1 K 



2 



dr] d 2 i] 



(24) 



Finally, we take into account the damping effects of the 
viscous forces acting on the antenna through the virtual 
work 
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2N 
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dx J d^ C v 



dy_ <V 

at at. 



S V , (25) 
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where C y and C v are the viscosities of the beam and the 
cantilevers, respectively. 

To account for the fact that damping effects and driv- 
ing forces are of the same order of magnitude as the non- 
linear effects, we scale the viscosities as C y = e 2 c y and 
C n = e 2 c v and the forces as Fj = e 3 fj, j = 1,2. We 
now proceed as explained in Ref. (l3l | to derive a time- 
averaged Lagrangian from Eqs. (f2"T|) , (j2"2"|) and (f2"I| . We 
substitute (f2Tj) into the Lagrangian (|22"|) and also into 
the virtual work (|25[) . perform the spatial integrations 
and keep the slowly varying terms only - i.e. those that 
are either constant or function of T 2 only. This yields: 



= {Mm (AjA'j - A'jAj) + C J3 \Aj 



+T 3 (Aj 



cc)}+2C 12 \A 1 \ 2 \A 2 \ 2 (26) 



(SW) 

where 
M 



QjSAj + cc, with Qj—2iujj^ijAj, (27) 

i=i 



3 - m b + 2Nm t Ljj ; T 3 = 77/^4 



-(Lc y + 2Nla n A j j) 



c _ T i™bUJ 2 



3T 2 £ b 



2AT, 



L 2 L 5 Z 2 

r!m b (w 2 6T 2 £ b 2NT. 



73 -"-jjjj 



L 2 



L 5 



Z 2 



-■yg(ifll22 + ^2211 F 4i^i2l2) 



and where the coefficients Tj, Lij, A^, I^j and -fQifcz are 



given in appendix VK\ 



mt(w\ + u>l)Ii 2 
(28) 

ijkl 



and separating the real and imaginary components of 
Eqs. PJJ), yields 



Wifa - flJ)Af i0i = j (C« a 2 + C12 a 2 ) - T cos(^), 
ojiiMia'i + HiOi) = J^sin(^i). (32) 

Looking for steady state (periodic) solutions, we impose 
a'j = and O'j = anc 
amplitude relations as 



a'j = and 9'j = and we finally obtain the frequency 



ai = iJ^^^+^^^VM^f - ^' (33) 



02 = 



together with 
tan6>i = 



[C 22 a*+C 12 a 2 ] ±^\ T p^-^. (34) 



Ca of + C12 a 2 - 4MjW i o- i 



i ^ j. (35) 



In first approximation the steady state solution can be 
cast into the form 



1 



1 X2(«) J COS ( W< ~ 6,2 )j • ( 36 ) 

The amplitudes a\ and a 2 are assumed small enough for 
the perturbation expansion to hold (notice that the book- 
keeping parameter e has been absorbed in the amplitudes 
and that Eqs. (|33|) - (p4)l can be used as such provided the 



dctunings <jj are redefined as <7j = s aj, the viscosities 
fij as = e 2 \ij and the forces Tj as Tj = s 3 T 3 ). 



D. Discussion 



C. Modulation equations 

Applying the extended Hamilton principle (see Ref. 
[HI]), we obtain the equations of motion for the modu- 
lations Ai(T 2 ) and A 2 (T 2 ) as 

From Eqs. and (f2"T)l , we then derive the following 

pair of modulation equations 



2iuj l (M i A' i + fnAi) = -2A % (C u \Ai\ 2 + 



+ C 12 \Aj\')+Tie^, (30) 

where G {1,2}, i ^ j. Looking for solutions in the 
polar form 

MT 2 ) = l - ai {T 2 )e l ^- e ^\ i e {1,2}, (31) 



The frequency- amplitude relations (|3"3"1) - (|3"4")) allow us 
to evaluate the frequency shift of the fundamental peak 
induced by a driving of the higher mode at the exact 
linear resonance frequency, lo 2 . This frequency shift is 
determined as the difference between the maximum of the 
resonance peak of the fundamental mode and u}\. Now, 
the amplitude a± becomes maximum if the square root in 
the r.h.s. of (f3"3"| vanishes, that is, 



T 



(37) 



For a system whose higher mode is driven exactly at fre- 
quency to — lo 2 we have of course a 2 — and then the 
amplitude of the second peak is solution to 



1 



AM 2 u> 2 



C22 a 2 + ^12 



2 2 
fifujf 



± 



F 2 



j4_ 

M 2 u\a\ M| 



= 0, 



(38) 
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which is a cubic equation for a 2 . Once the solution 
a 2 {T\, T 2 ) is known, we can reinstate it in (|33[) and we 
obtain the frequency shift, 01 = Q — lui as 



a 1 (Ti , ^2 



1 



4MiWi 



jp2 

C11 — 2 — 2 



Cl2 a 2 (-^1) ^2) 



which provides the frequency shift as a function of the 
forces (or driving power), Ti and T 2 , of the fundamental 
and excited modes. 



V. CONCLUSION 

Here we have expanded and generalized the analyti- 
cal modeling of coupled-resonator elastic structures. Our 
one-dimensional continuum model is shown to capture 
the relevant dynamics and response spectrum of antenna- 
like structures with arbitrary (symmetric) cantilever dis- 
tributions as well as boundary conditions on the central 
beam. A detailed analysis of the band structure of the 
response spectrum in the special cases of constant and 
alternating cantilever length arrays resonator provided 
insights into the modal dynamics of the cantilevers ver- 
sus the carrier beam. Design parameters such as can- 
tilever lengths are shown to modify the response spec- 
trum and should help achieve potentially advantageous 
performance. 

We further showed that perturbation theory can be 
used to analyze the nonlinear modal response of a pro- 
totype antenna resonator. The analysis was performed 
without loss of generality on two modes that have been 
selected for their practical advantages in performance. 
The anharmonic amplitude-frequency coupling that re- 
sults from the inclusion of nonlinear terms in the equa- 
tions of motion elucidates the experimental behavior of 
hybrid nano-resonators under strong excitation. Thus the 
current analysis provides a solid analytical foundation to 
the understanding of device performance, with direct ap- 
plication to a wide range of current resonator devices, as 
well as a forthcoming experimental report by the authors. 
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APPENDIX A: EFFECTIVE COEFFICIENTS 

We provide here the definition and the numeri- 
cal/analytical evaluation of the quantities involved in 
Eqs. (j2"5|) for the antenna dimensions given in 0]. 



r x 



[y>i(f)] dv\ du, 



\jq 
1 



<Pl(u) du, 







r± = I fi(u)du, 
Jo 

Lij = / hi(v)hj(v) dv, 



A 







hi(v)(hj(v) - 1) dv, 



lij — 



^ (J^h^h^dv) dv, 

K ; ,u ^ / h'iWtijWhttvWWdv. (A2) 




r 2 = / [<Pi{u)<p"{u)] du, 



(Al) 



with hj(v) = Xj( v ) + !• All the integrals above can nor- 
mally be evaluated analytically. For Ti, i £ {1, ... ,4}, 
results are simple enough to be displayed below. 



Ti = ^{0it + 2) ~ 6.1513, 

T 2 = ^(5/3ii+ 11) - 2846.4975, 

At 

r 4 = — ~ -0.8308 

Pi 



where t = tan(/?i/2). Note that cos(/3i) = 1/ cosh(/?i) 
and sin(/?i) = — tanh(/3i). Apart from A;,; = La — 
2^2(7i)/7i and La = £(71^), the other integrals have 
been evaluated numerically for the parameters given in 
7] . We have found 

L n ~ 1.00012 ; L 22 ~ 3.89887, 

A n ~ 6.16 x 10~ 5 ; A 22 ~ 4.9687, 

I n ~ 1.64 x 10~ 16 ; 7 22 ^ 232.49, 

I12 ~ 1.94 x 10~ 7 ; K 1212 ~ 8.26 x 10~ 7 , 

Kim - 6.58 x 10~ 16 ; ^2222 ^ 1-07 x 10 3 , 

Ku2 2 ~ 9.84 x 10~ 7 ; K 2211 ~ 6.98 x 10~ 7 . 

So that, finally, 

Mi = 1.74 x 10" 14 ; M 2 = 4.17 x 10 -14 , 
Ti = -4.44 x 10~ 6 /i ; C l2 = 1.65 x 10 17 , 
Cn = -5.07 x 10 13 ; C 22 = 1-05 x 10 21 , 
uji = 1.55 x 10 s ; uj 2 = 1.85 x 10 10 . 

where all the results are given in SI units. From the quan- 
tities ujj, the frequencies of the fundamental and collec- 
tive modes are 24.7 MHz and 2.94 GHz, respectively. 
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